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ABSTRACT 

We use the Millennium Simulation, a large, high resolution N-body simulation of the 
evolution of structure in a ACDM cosmology, to study the properties and fate of 
substructures within a large sample of dark matter haloes. We find that the subhalo 
mass function departs significantly from a power law at the high mass end. We also find 
that the radial and angular distributions of substructures depend on subhalo mass. 
In particular, high mass subhaloes tend to be less radially concentrated and to have 
angular distributions closer to the direction perpendicular to the spin of the host halo 
than their less massive counterparts. We find that mergers between subhaloes occur. 
These tend to be between substructures that were already dynamically associated 
before accretion into the main halo. For subhaloes larger than 0.001 times the mass 
of the host halo, it is more likely that the subhalo will merge with the central or 
main subhalo than with another subhalo larger than itself. For lower masses, subhalo- 
subhalo mergers become equally likely to mergers with the main subhalo. Our results 
have implications for the variation of galaxy properties with environment and for the 
treatment of mergers in galaxy formation models. 
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1 INTRODUCTION 

The presence of substructures within dark matter haloes is a 
distinctive signature of a universe where structures grow hi- 
erarchically. Low mass objects collapse at high redshift, and 
then increase their mass by smooth accretion of dark mat- 
ter or by merging with other haloes. Once a halo is accreted 
by a larger one, its diffuse outer layers are rapidly stripped 
off by tidal forces. However, the core, which is much denser, 
generally survives the accretion event and can still be rec- 
ognized as a self bound structure or subhalo within the host 
halo for some period of time afterwards. 

In early N-bo dy simulations, ha lo es ap peared as fairly 
smooth objects (|Frenk et al.l Il985l . Il988t ). However, as 
the attainable mass and force resolution has increased, 
subhaloes have been identified and their properties studied 
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population have important implications for galaxy forma- 
tion, dark matter detection experiments and weak lensing. 
For instance, subhaloes are expected to host satellite 
galaxies within groups and clusters and their evolution 
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once inside the host could give rise to observable changes. 
In particular, a merger between two substructures could 
trigger an episode of star formation or a mo rphological 
transformation (e.g. ISomerville fc Primac^ll999t ). 

In spite of this, the merger history of subhaloes remains 
relatively unexplored. This is a challenging problem which 
demands a simulation with high mass and force resolution. 
In particular, obtaining a statistical sample of mergers in- 
volving the largest substructures requires a large sample of 
host haloes. Most studies of substructure in halos have fo- 
cused on resimulating, at very high resolution, a small num- 
ber of halos selected from a larger, lower resolution simu- 
lation. However, by studying only a few haloes, important 
aspects related to variations produced by differences in the 
accretion and merger histories of haloes, as well as any in- 
fluence of the environment, could remain hidden. This ap- 
proach may also introduce systematic biases arising from the 
criteria used to select the haloes to be resimulated. 

In this paper, we overcome these problems by using the 
largest dark matter sim ulation published to date, the Mil- 
lennium Simulation (MS. ISpringel et alllioosh . The MS pro- 
vides a large cosmological sample of dark matter haloes and 
associated substructures spanning a considerable range in 
mass, allowing us to assess robustly the properties and fate 
of the subhalo population. We complement our results with 
a higher resolution simulation of a smaller volume (hereafter 
HS) which has a particle mass almost ten times smaller than 
that used in the MS (Jenkins et.al, in prep). 

The layout of this paper is as follows. In Section (2] we 
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briefly describe the simulations used in this work along with 
the properties of our halo and subhalo catalogues. In Sec- 
tion |3] we investigate some general properties of subhaloes, 
namely their mass function, radial distribution and spatial 
orientation with respect to their host halo. The exploration 
of substructure mergers and destruction is presented in Sec- 
tion |4l Finally, we summarize our findings in Section [S] 



2 METHOD 

In this section we describe the N-body simulations we have 
analyzed in this work. We also discuss the identification and 
characterization of the halo and subhalo catalogues. 



2.1 N-body Simulations 

The main simulation on which our analysis is based is the 
Millennium Simulation (|Springel et al.ll2005l ). The MS cov- 
ers a comoving volume of 0.125 /i~^Gpc^ of a ACDM uni- 
verse in which the dark matter component is represented 
by 2160^ particles. The assumed cosmological parameters 
are in broad agreement with those der ived from joint anal - 



yses of the 2dFGRS galaxy clustering (IPercival et al.l 2001 



and WMAPl microwave background data ( Spergel et al.l 
l2003l : ISanchez et alll2006l). as well as with those derived 
from WMAP5 data (|Komatsu et all 120081 ). In particular, 
the total mass-energy density, in units of the critical den- 
sity, is f2ni = f^dm + f2b = 0.25, where the two terms 
refer to dark matter and baryons, with S7b = 0.045; the 
amplitude of the linear density fluctuations in 8/i~^Mpc 
spheres is erg = 0.9; and the Hubble constant is set to 
h = ffo/(100kms"^Mpc"^) = 0.73. The particle mass is 
m-p = 8.6 X 10* h~^MQ and the Plummer-equivalent soften- 
ing of the gravitational force is e = 5/i~^kpc. 

To complement our results and to test for numerical 
effects we have also employed another simulation with bet- 
ter mass resolution to which we refer as HS. This simula- 
tion follows 900^ dark matter particles in a ACDM cube of 
side 100/i~^Mpc. The HS assumes the same cosmological 
parameters as the MS. However, the smaller box yields a 
smaller particle mass, nip = 9.5 x 10^ Mq, so objects of 
a given mass are resolved with almost 10 times more parti- 
cles than in the MS. Finally, in the HS the softening length 
is e = 2.4ft,~^kpc. 

The MS and HS were carried o ut using a me mory effi- 
cient version of the Gadget-2 code (|Springelll2005l '). 



2.2 Halo and Subhalo catalogues 

In both simulations, particle positions and velocities are 
written at 64 output times which, for z < 2, are roughly 
equally spaced in time by 300 Myr. In each of these outputs 
we have identified dark matter haloes usin g the friends-of- 
friends (FoF) algorithm (|Davis et al.lll985l ). with a linking 
length of 0.2 times the mean interparticle separation. The 
volume and particle number of the MS provide a unique 
resource of well resolved haloes to study. By way of illustra- 
tion, there are 90891 haloes ai z — with mass in excess 
of 5.4 X 10^^ Mq (one of the bins we use below), which 
corresponds to 6272 particles; ai z = 1 the number of haloes 
in excess of this mass is still 61481. On the cluster-mass 



scale, for example, there are 356 haloes at z = which are 
more massive than 4 x 10^* Mq, corresponding to 464576 
particles. 

Well resolved FoF haloes are not smooth, but contain 
a considerable amount of mass in the form of substructures. 
These dark matter clumps or "subhaloes" are identified and 
catalogued using a modified version of the subhalo finder 
algorithm, SUBFIND, originally presented in ISpringel et al.l 
(,2001 ). The statistics of the subhalo catalogue are impres- 
sive. At z = SUBFIND fists 339840 structures with more 
than 200 particles in the MS within haloes of at least 
5.4 X 10^^ H-'^Mq. M z = l there are 194629 substructures 
with the same characteristics. Note that SUBFIND not only 
finds substructures within a FoF halo, but it is also capable 
of identifying substructures within substructures. 

An important issue for studies of substructures is the 
definition of the boundary and position of the host halo. In 
our analysis, the centre of the halo is defined as the position 
of the most bound particle (i.e. usually the one possessing 
the minimum gravitational potential). This choice for the 
halo position agrees, within the softenin g length, with that 
found by a shrinking sphere algorithm (|Power et al.ll2003h 
for 93% of the haloes that are resolved with 450 or more 
particles. As shown bv iNeto et al.l (|2007l ). the 7% of cases 
in which the two methods disagree are due to the FoF algo- 
rithm artificially linking multiple structures. In these cases 
the position of the most bound particle provi des a more 
robust identification of the centre, as noted by iNeto et al.l 
(|2007l ). 

We define the halo boundary as the sphere of radius 
r2oo which contains a mean density of 200 times the critical 
density, pcrit. Therefore, the mass of the halo is: 
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We keep in our catalogues only subhaloes within this 
sphere. Although the choice of the factor of 200 is motivated 
by the spherical collapse model in a Einstein-de-Sitter uni- 
verse, it is somewhat arbitrary for our ACDM simulations. 
However, the r2oo definition has the advantage of being in- 
dependent of both redshift and cosmology. Moreover, it has 
became a de facto standard in studies of substructures. Nev- 
ertheless, we have tested our results against other definitions 
of the halo boundary without finding any qualitative differ- 
ences. In the following, when we refer to the mass and radius 
of a host halo, we always mean M200 and r2oo. 

Finally, we build m erger trees using an algorithm simi- 
lar to that described bv lSpringel et akl l|2009' ) which follows 
the evolution of subhaloes. In this way, we can identify the 
haloes and subhaloes that will be involved in a merger during 
a subsequent snapshot. Note that these merger trees are con- 
structed using only the information contained in the FoF and 
SUBFIND catalogues, and there is no attempt to force mass 
conservation, as would be required if the merger trees were t o 
be used in a galaxy formation code fsee lHarker et al1l2006l ). 
The descendant of a subhalo is defined as the structure that 
contains the majority of the 10 percent most bound parti- 
cles from the subhalo. When two satellite subhaloes have 
the same descendant in a following snapshot, we tag such 
an event as a substructure merger. 
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Figure 1. Top row. Differential number of substructures per host halo as a function of their mass relative to that of the host halo, 
A^aub/^host- Note this is the mass of the subhalo at the redshift labelled, in some cases after substantial stripping of mass has taken 
place. Solid lines show the results from the MS while dashed lines show the results from the HS. In both cases, lines of different colours 
show the subhalo mass function in host haloes of different masses (as indicated by the legend). Each column shows a different redshift, 
as labelled. At each redshift, the dotted lines display the overall best fit of our model, Eq. [2] with the parameter values given in the 
legend. Parameters of the fits to individual mass bins at z = are listed in Table 1. Bottom row: Relative difference between the overall 
best fit and measurements of the subhalo mass function for different host masses. The dot-dashed line shows the difference between our 
model, Eq. (2] and a power-law fit. Only results for subhaloes which are resolved with more than 50 particles are shown. 



3 SUBHALO PROPERTIES 

Before presenting our results regarding subhalo mergers, we 
consider some general properties of the subhalo population. 
Although some of these properties have been studied by pre- 
vious authors, the large volume and high resolution of the 
MS and HS reveal some features which were inaccessible 
to earlier work. Furthermore, the knowledge of the subhalo 
properties will help us to understand the results presented 
in the next section. 



3.1 Subhalo mass function 

We first consider the distribution of subhalo masses - the 
subhalo mass function. The top panels of Fig. [T] show the 
mean number of substructures within dark matter haloes, 
per host halo, per logarithmic interval in subhalo mass. The 
results are displayed as a function of subhalo mass relative 
to the mass of the halo in which it resides, Msub/A/host. In 
this way we can easily compare results across a range of 
halo masses. In the ranges of overlap, the results from the 
MS and HS agree well; this provides a useful, but limited, 
test of the convergence of our results. 

For the redshifts plotted in Fig. [l] there is only a small 
variation of the subhalo mass function with host halo mass. 
Indeed, a universal function describes the behaviour reason- 
ably well over the range of subhalo mass resolved by our 
simulations: 



dN 



dln(Msub/Mho 



A 



exp 



Afhost 
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where A'^ is the number of subhaloes per host halo. The val- 
ues of A, a and a in this overall fit at each redshift are given 
in the legend of Fig. [1] For this overall fit, we have forced 
the slope a to have the same value independently of redshift. 
In general, we find that a — —0.9 is a good approximation 
to the best fit from z — to z = 2.5. It is also important 
to n ote that the powe r-law fit widely used in the literature, 
(e.g. iGao et aLll2004 ) is only valid over a limited range of 
fractional subhalo masses, Msub/Mhost < 0.04. We also see 
that the maximum subhalo mass for which a power-law suc- 
cessfully describes the mass function decreases at higher red- 
shifts, Afsub/Mhost ~ 0.015 at 2 = 1 and Msub/Mhost ~ 0.04 
at 2 = 0. The bottom panels of Fig. [T] show the relative 
difference between the fit given by Eq. ((2)| and the mass 
function of subhaloes measured in host haloes of different 
masses. 

We have also fitted Eq. [5] to the subhalo mass func- 
tions in each halo mass bin, this time letting the slope a 
vary; we list the best-fit parameters for jz = in Table [1] At 
the low fractional mass end, where the subhalo mass func- 
tion behaves as a power-law, we generally find slopes that 
are lower than the critical value, a = —1 (which separates 
divergent from convergent mass functions). The slopes we 
find are in broad agreement with previous estimates of the 
power-law index of the subhalo mass function, which range 
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a 






9.2 X 1012 


-2.05 


-0.87 


0.17 


-1.8 




2.7 X 10^^ 


-2.06 


-0.89 


0.16 


-2.5 


MS 


7.9 X 10^^ 


-1.98 


-0.88 


0.13 


-2.8 




2.3 X lO^-* 


-2.00 


-0.90 


0.10 


-3.5 




6.8 X lO^-* 


-1.86 


-0.87 


0.06 


-3.8 


HS 


3.1 X 10^2 

9.2 X 10^2 
2.7 X 10^3 


-2.00 
-2.05 
-2.15 


-0.83 
-0.88 
-0.93 


0.16 
0.17 
0.14 


-2.5 
-2.8 

-3.5 



Table 1. The bost-fit parameters to the mass function of sub- 
haloes residing in haloes of different mass at z = 0, using Eq. [2] 
The columns are as follows: (1) The N-body simulation from 
which the halo sample was extracted. (2) The mean mass of the 
host haloes. (3) The logarithm of the amplitude. (4) The power- 
law index. (5) The damping strength. (6) The minimum fractional 
subhalo mass included in the fitting. 



from -0.8 to -1.0 (iMoore et a.l.lll999l: iGhigna et al-llioOC 



De Lucia et"aIll2004inG'ao et al.ll2004l: IPiemand et al.ll2004 : 
Shaw et al .l2007l :l iDiemand et al .l2007f ). In particular, our re- 
sults agree with those from the much higher resolution sim - 
ulations of individual galactic halos of lSpringel et al.l (|2008l l , 
but are inconsistent with the steeper slope advocated, also 
for galactic halos, bv lDiemand et ahl (|2008| ). 

At the high mass end, the subhalo mass function de- 
parts from a power-law and decreases exponentially. This 
behaviour was previously detected in N - body simulations 



(at lower significance ) by Giocoli et al.l (ioosi) (and pre- 
dicted analytically bv lvan den Bosch et al. I I2OO5I I. However, 



this feature was not apparent in earlier studies which used 
resimulations of individual haloes. Resimulations of single 
objects have the advantage that computational effort can 
be focused. A halo can be resolved with a vast number of 
particles and its substructures identified over a large range 
of masses. Unfortunately, this approach comes at the price 
of losing the rich information contained in the variety of as- 
sembly histories, relaxation states and, more importantly, 
the population of high mass subhaloes. As can be seen from 
Fig. [1] the abundance of these objects is much lower than 
that of smaller subhaloes - usually we would find just a few 
in each halo. Because these halos are so rare, the damping 
of the power- law at high i\fsub/Afhost is missed in individual 
resimulations. By contrast, with the huge sample of haloes 
and their massive subhaloes in our analysis, we can robustly 
probe this subhalo mass range. 

Even tho ugh the subhalo ma ss function appears roughly 
universal (e.g. iMoore et al.1119991 1 , we have detected at every 
redshift a small dependence on the mass of the host system. 
Small substructures of the same fractional mass are more 
abundant in high mass haloes than in low mass haloes. This 
correlati on has also been seen in a number of other stud- 
ies (e .g. iGao et all |2004 IShaw et all l2007l : iDiemand et~all 
I2OO7I ). However, we also find evidence that this trend holds 
only in the power-law region of the subalo mass function and 
actually reverses at the high mass end - low mass haloes seem 
to host relatively more massive subhaloes than do high mass 
haloes. 

Perhaps surprisingly, the variety of features present in 



the mass function of subhaloes is consistent with a relatively 
simple picture. There are two key ingredients that shape the 
subhalo mass function: (i) the mass function of infalling ob- 
jects and (ii) the dynamical evolution of subhaloes orbiting 
within the host halo due to dynamical friction and tidal 
stripping. The first of these is responsible for the univer- 
sality described above and sets t he subhalo ma s s func tion 
to first order. As first found by iLacev fc Cold (|l993l ) us- 
in g the extended Press Schechter formalism, and confirmed 
bv lGiocoli erai] l|2008l ) using N-body simulations, the mass 
function of subhaloes at infall is almost independent of host 
halo mass and redshift when expressed as a function of 
Msub/Mhost, and can be described as a power-law with a 
high mass cut-off. 

After subhaloes fall into a host halo, their orbits sink 
due to dynamical friction and, at the same time, the sub- 
haloes lose mass due to tidal stripping. These processes 
cause the subhalo mass function to evolve away from its 
form at infall. The rates for these processes depend on the 
fractional mass of the subhalo, Msub /Mhost , and on the dy- 
namical timescale of the host halo. Therefore, if all haloes 
had identical structure and assembly histories, these pro- 
cesses would preserve a universal form for the subhalo mass 
function, independently of Mhost • However, haloes of differ- 
ent masses on average assemble at different redshifts in spite 
of the similar mass function of subhaloes at infall, and this 
breaks the u niversal shape of the subhal o mas s function, as 
discussed by Ivan den Bosch et al.l (|2005l ) and iGiocoli et al.l 
(j2008h . On average, massive haloes are younger than their 
less massive counterparts a nd they are more li kely to have 
experienced recent mergers (|Lacev fc Colelll993l ). These pro- 
vide a fresh source of substructures which have had less 
time for orbit decay due to dynamical friction and to be 
tidally stripped. High mass haloes are therefore expected to 
have more substructures than low mass haloes. Another ef- 
fect which acts in the same direction is that small haloes 
tend to accrete their subhaloes at higher redshifts when dy- 
namical timescales are shorter. As a result, they strip out 
mass from the substructures more quickly than large haloes, 
where massive substructures can survive for longer. 

3.2 Most massive subhaloes 

The high-mass tail of the distribution of substructure is 
examined in greater detail in Fig. [5] The three panels in 
this plot display the distribution of the fractional mass, 
Afsub/A/host , for the first, second and third largest substruc- 
tures within haloes of different mass at 2; = 0. As before, 
results from the MS and HS agree very well. 

In contrast to the results presented in the previous sub- 
section, the distributions of fractional masses seem to be 
independent of the host halo mass. (We have also checked 
that they are independent of redshift.) In particular, in every 
halo, the fractional masses follow a log-normal distribution 
with mean (log^o (Msub/Mhost)) = -1.42, -1.79 and -1.99, 
and standard deviation o'iogio(A4ub/Mhost) ~ 0.517, 0.382 and 
0.348 for the each of the three largest subhaloes respectively. 
Albeit with considerable scatter, these values imply that the 
most massive substructure contains typically 3.7% of the to- 
tal mass of the halo while the second and third most massive 
subhaloes contain 1.6% and 1% of the mass respectively. 

Due to the large dispersions, the distributions can only 
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Figure 2. The distribution of the fractional mass, Msub/Mjjost, 
of the 1st, 2nd and 3rd largest substructures in haloes of different 
mass at 2 = 0. The solid Hues show the resuhs from the MS vfhile 
the dashed lines show the results from the HS. In each panel, the 
black solid lines indicate the log-normal function that best fits our 
results. Note that only substructures resolved with 20 particles or 
more are displayed. 



be measured reliably in haloes resolved with a large number 
of particles. For instance, the mean fractional mass of sub- 
haloes is overestimated for haloes resolved with fewer than 
~ 1000 particles (the exact limit depends on the scatter and 
mean of the true distribution) , i.e. ~ 1 X 10^^ h-^MQ in the 
MS and ~ 1 x lO" h~'^ Mq in the HS. The upward bias is 
caused by the finite resolution of the simulations (there is 
a limit on the smallest subhalo that we can identify) which 
truncates the low mass tail of the distribution of fractional 
masses. 

Hints of a universal behaviour of the fractional 
masses of the l argest subhaloes were already detected by 
Lucia et al.l l|2004 ) (although they claim a weak depen- 
dence with host halo mass). Our results are broadly con- 
sistent with theirs but, with the large halo catalogues from 
the MS and HS, we are able to probe the full probability 
distribution function robustly. 

The apparently universal shape of these distributions 
could, in principle, be understood within the broad pic- 
ture just discussed. Presumably it reflects the distribution 
of masses of the infalling haloes whic h, as we have seen , 
is independent of the host halo mass l|Lacev fc Cold 1 19931 : 



iGiocoli et al.|[2008l ) . The large scatter must then result from 
the large range of accretion histories at a given host halo 
mass. We leave further investigation of these ideas for fu- 
ture work. 

3.3 Radial distribution of subhaloes 

Fig.|3]shows the number density of subhaloes as a function of 
radius, relative to the mean number density of substructures 
within r2oo in the same fractional mass range. Each panel 
focuses on substructures of different masses, from small sub- 
haloes (10"" < Msub/Mhost < 10~^) in the leftmost panel to 
large ones (10~^ < Mgub/Afhost < 1) in the rightmost panel. 
As in previous plots, lines of different colours show results 
for subhaloes that reside in haloes of different mass, and the 
different line types (solid and dashed) indicate the results 
for the two simulations. We also plot the radial profile of 
the dark matter as a black dotted line in each panel. 

Comparison of the MS and HS indicates that our re- 
sults are insensitive to the mass resolution (although the 
overlap between th e two simulation s is limited). As in pre- 
vious studies (e.g. iGao et ahl |2004| ). we find that the ra- 
dial distribution has little dependence on the host halo 
mass at a given A/sub /Mhost- This is quite remarkable since 
each panel mixes subhaloes that: (i) are resolved by num- 
bers of particles that differ by orders of magnitude and (ii) 
occupy haloes which are in a variety of dynamical states 
(age, relaxation, etc). We also see that in all cases, the 
radial distribution of subhaloes is less centrally concen- 
trated than the dark matter, as was also found in previ- 



ous studies (e. 
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In addition, we see a significant difference between the 
distribution of massive subhaloes (Msub > 10~^Mhost) and 
that of small ones {Msuh < 10^'^Mhost)- While the overall 
radial profiles seem to be fairly independent of subhalo mass, 
the more massive subhalos tend to avoid the central regions 
of the host halo, while the less massive on es have a more 
centr ally concentrated distribution (see also lDe Lucia et al.l 
l2004h . H owever, the distributi ons agree in the outer parts of 
the halo. ISpringel et alJ l|2008l ) found a similar effect to ours 
in the Aquarius set of simulations of galactic halos which, 
although limited in number, span a huge dynamic range in 
subhalo mass. 

These dissimilar density profiles for different sub- 
halo masses have a simple dyna mical explanation (e.g. 
iTormen et "all 1 19981 : iNagai fc Krav tsov 2005). Once a halo 
falls into a more massive system, dynamical friction and 
tidal striping start to act. The accreted subhalo will rapidly 
be stripped of its outer layers and will lose a significant 
fraction of its mass during the first pericentric passage. This 
mechanism naturally differentiates the radial distribution of 
substructures of different masses: massive structures sink 
more rapidly due to dynamical friction and, as a result, also 
lose mass more quickly by tidal stripping. Therefore they do 
not survive long in the central regions, in contrast to small 
subhaloes. The massive subhaloes which are present in the 
halo must have been accreted more recently than the aver- 
age low mass subhalo. The timescale for dynamical friction 
depends on the relative mass of the subhalo and its host 
halo, not on their absolute values, which would explain the 
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Figure 3. The number density of subhaloes relative to the mean within -Rhosti a function of the distance to the centre of their host 
halo, in units of the radius of the host halo, -Rhost • Each column corresponds to a different fractional mass range for subhaloes, while the 
rows display the results at two separate redshifts. The solid and dashed lines show the density profiles from subhaloes in the MS and HS 
respectively. In both cases, the different colours correspond to subhaloes residing in haloes of different mass as shown in the legend. The 
black dotted lines in each panel indicate the mean dark matter density profile of haloes in our simulations. Results are shown only for 
subhaloes resolved with at least 200 particles. The top row shows results for z = 0.5 and the bottom row for 2 = 0. 



approximate independence of the distribution on the host 
halo mass. 



3.4 Angular distribution of subhaloes 

To end this section we investigate the angular distribution of 
subhaloes within dark matter haloes. Previous work has ex- 
amined the relationship between the angular distribution of 
substructures and the shape of the host ha lo (|Tormenll 19971 : 
iLibeskind et al ] l2007l : iKnebe et al.l l2008al lbh. Here, we ex- 
amine instead the orientation relative to the spin axis of the 
host halo. Fig.Ushows the probability distribution function 
of the cosine of the angle between the angular momentum 
vector of the host halo and the vector joining its centre with 
that of the subhalo. We show results for two separate ranges 
of subhalo mass: subhaloes with mass smaller than 2% of 
the host halo mass (dashed lines) and those with masses 
greater than 2% (solid lines). We distinguish different host 
halo masses by different colours, and show different redshifts 
in different panels. Note that we only display results for the 
MS simulation for clarity. 

As shown by Bett et al. (2007), the accuracy of the mea- 
surement of spin direction in the MS degrades significantly 
(uncertainty > 15deg) for haloes resolved with fewer than 



1000 particles or for those where the spin magnitude, is 
such that: 

< lo-l■^ (3) 

V G A/host -Rhost 

where G is Newton's gravitational constant. Although the 
inclusion of haloes that do not satisfy these criteria does not 
seem to affect our results quantitatively, we have chosen to 
show only those haloes that met these requirements, so that 
the angle relative to the spin axis can be reliably determined. 

We see from Fig. |4] that the angular distribution of 
subhaloes tends to be aligned perpendicular to the spin 
axis of the host halo. (We remind the reader that in this 
plot, an isotropic angular distribution would correspond to 
a horizontal line, while a distribution aligned at 90deg to 
the spin axis will peak around cos 6 ~ 0.) The strength 
of this alignment effect depends on the fractional subhalo 
mass, Maub/Mhoat, being much stronger for higher mass sub- 
haloes. We also see that the angular distribution for a given 
Msub/Afhost is almost independent of the host halo mass and 
the redshift (see also[Kang ct al. 2001). 

We can understand this behaviour qualitatively as re- 
flecting the growth of haloes by the accretion of dark matter 
(in halos or more diffuse form) along filaments. The central 
regions of haloes acquire most of their angular momentum 
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Figure 4. The probability density distribution of the cosine of 
the angle 8 between the angular momentum vector of the host 
halo and the vector joining its centre with that of the subhalo. 
Each panel shows results for a different redshift, as indicated by 
the label. The lines in each panel display the distribution for 
subhaloes in two different mass bins: 0.02 < Mgub/Mjjost < 1 
(solid lines) and 0.004 < Afgub/Mhost < 0.02 (dashed lines). Linos 
of different colour indicate subhaloes residing in host haloes of 
different masses, as shown in the legend. An isotropic angular 
distribution corresponds to a horizontal line. 

at a relatively late stage from the orbital angular momen- 
tum of this infalling material, and so will tend to have spin 
axes perpendicular to the curr ent filament (e.g IShaw et al.l 
l2006l : lArag6n-Calvo et ai]l2007l ) . On the other hand, insofar 
as the subhaloes "remember" the direction from which they 
fell in once they are orbiting inside the host halo, then their 
spatial distribution will tend to be aligned with the filament 
from which they were accreted, and so will be perpendicular 
to the spin axis. We can also understand the dependence 
of the strength of this alignment on subhalo mass in this 
picture. Subhaloes with large Msub/Mhost on average have 
been orbiting in the host halo for less time than haloes of 
lower Msub/Mhost, due to the combined effects of dynam- 
ical friction (which causes higher mass subhaloes to sink 



faster) and tidal stripping (which converts high-mass sub- 
haloes to low mass). We expect subhaloes increasingly to 
lose memory of their initial infall direction the longer they 
have orbited in the host halo (which in general is lumpy and 
triaxial). Since high Msub/Mhost subhaloes have undergone 
fewer orbits, their current angular distribution should be 
more closely aligned with their infall direction, and there- 
fore with the current filament, compared to subhaloes of 
lower mass. 

Our results seem generally consistent with previous sim- 
ulation results on the alignment of the subhalo distribution 
with the shape of the host halo, and the relationship be - 
tween the shapes and the spin axes of halos. iTormenI (|l997l ) 
found that the angular distribution of subhaloes as they fall 
into a host halo (crossing through r2oo) is anisotropic, and 
tends to be aligned along the ma jor axis of the host halo . 
Previous studies (e.g.jKnebe et an i2004: Zo ntner et al.|[2005l : 
iLibeskind et al.ll2007l l found that the angular distribution 
of subhaloes within a host halo is aligne d along the major 
axis of the host halo. On the other hand. lBett et al.) (|2007l ) 
showed that the angular momentum of a halo is generally 
aligned with its minor axis and perpendicular to its ma- 
jor axis. Putting these results together, we would expect 
the subhalo distribution to be aligned perpendicular to the 
spin axis of the host halo, but ours is the first study to 
demonstrate this directly, and also to demonstrate that the 
strength of the alignment depends on subhalo mass. 



4 MERGERS BETWEEN SUBHALOES 

As we have seen, once a halo is accreted by a larger one, its 
outer layers are rapidly stripped by tidal forces. However, 
the core generally survives the accretion event and can still 
be recognized as a substructure or satellite subhalo within 
the host halo for some time afterwards. Furthermore, not 
only may the main infalling halo survive, but also substruc- 
tures within it. In this case, there are substructures inside 
substructures. 

While orbiting inside the halo, dynamical friction causes 
the orbit of a subhalo to lose energy and to sink towards the 
centre of the host halo. As the subhalo sinks, it suffers fur- 
ther tidal stripping. Eventually, the subhalo may be totally 
disrupted: there is a merger between the satellite subhalo 
and the central subhalo. Nevertheless, on its way to destruc- 
tion, a subhalo can survive for several orbits during which 
it may experience a merger with another satellite subhalo. 
In the following subsections we will investigate the merging 
of these substructures. 

The interaction between subhaloes was previously in- 
vestig ated in cosmological simulations by iTormen et al.l 
() 19981 ) ■ who studied the rate of penetrating encounters 
between satell i te su bhaloes, but not the merger rate. 
iMakino fc HutI ()l997l ) derived an expression for the merger 
rate between subhaloes in galaxy clusters based on an en- 
tirely different approach, motivated by the kinetic theory 
of gases. In this case, the merger rate per unit volume be- 
tween halos of mass Mi and M2 is i?mcrgo = ?^l^^2cr(^;l2)^'l2, 
where ni and n2 are the respective number densities, V12 is 
the relative velocity, and a-{vi2) is the merger cross-section. 
They used N-body simulations of isolated spherical halos to 
derive merger cross-sections for equal-mass halos as a tunc- 
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Figure 5. The mean number of satellite mergers per subhalo and per unit of time relative to the age of the universe, as a function of 
the mass of the progenitor of the less massive object involved in the merger. Two cases are displayed: the number of satellites destroyed 
or merging with the main substructure (top thin lines) and the number of mergers between two satellites (thick bottom lines). The solid 
lines show results from the MS while the dashed lines show results from the HS. As indicated by the legend, in both cases, coloured lines 
represent results for haloes of different mass. The three panels are for three different redshifts: z = 1.0, z = 0.5 and 2: = 0. Note that in 
the case of the merger between a satellite and a central structure, we show examples involving subhaloes of at least 200 particles, but 
we reduce the limit to 50 particles in the case of mergers between two satellites. In each panel, the legend states the redshift, the age of 
the universe, tn, and the time interval, dt, over which we measure the rates. 



tion of their relative velocity, and then assumed that merg- 
ers in clusters occurred between pairs of subhaloes drawn 
from random uncorrelated orbits, with a Maxwell ian dis- 
tribution of relative velocities. The lMakino"fc expres- 
sion was then extrapolated to the case of unequal sub- 
halo masses and incorpo rated into a semi-analytical m odel 
of galajcy formation by ISomerville fc Primack (|l999l ) and 
Hatton et aP (l2003h. We will investigate below whether the 
Makino &: Hud (|l997t ) kinetic theory approach has any ap- 
plicability to subhalo mergers in a realistic cosmological con- 
text. 



4.1 Subhalo merger rate 

Fig. [5] shows the mean merger rate of satellite subhaloes, 
plotted against the fractional mass of its progenitor. This is 
the mass of the satellite before accretion divided by the mass 
of the host halo at the time of the merger. The rate is nor- 
malized per subhalo, with time in units of the age of the uni- 
verse at that redshift. This normalized rate is thus roughly 
equal to the probability that a satellite subhalo will merge 
over one Hubble expansion time. A rate higher than unity in- 
dicates that the process happens on a timescale shorter than 
a Hubble time. There are two sets of curves in this figure: (i) 
the thinner, higher amplitude lines which show mergers be- 
tween a satellite and a central subhalo, as a function of the 
subhalo mass, and (ii) the thick lines which correspond to 
satellite-satellite mergers, plotted as a function of the mass 
of the smaller subhalo. As in previous plots, difi'erent line 



colours show different host halo masses, and different line 
styles (solid and dashed) show the two simulations used. 

We see from Fig. [S] that over most of the subhalo mass 
range resolved by our simulations (for Msub/Mhost 10"'^), 
it is more likely for a satellite subhalo to merge with the 
central subhalo than with another more massive satellite 
subhalo. For instance, at z = 1, taking into account all host 
haloes, there are 17155 satellites which merge with a central 
subhalo over one timestep, while the number of satellites 
involved in a merger with another satellite over the same 
period is 509, a ratio of 40 : 1. The situation is similar at 
z = even though the ratio decreases to 6 : 1 (1645 vs 290). 
In general, the likelihood of both merger rates slightly de- 
creases at lower redshifts. This may reflect the slower build- 
up of structure (relative to the Hubble time) as the universe 
becomes dominated by vacuum energy. 

As we consider smaller subhalo masses, we see a de- 
crease in the destruction rate (see the appendix for a dis- 
cussion of overmerging eff'ects due to insufficient mass res- 
olution). This may be due to the inefficiency of dynam- 
ical friction for low mass structures. On the other hand, 
there is an increase in the satellite-satellite merger rate as 
the subhalo mass decreases. Presumably this is due to the 
increasing number of potential merger partners, reflecting 
the form of the subhalo mass function. Additionally, the 
abundance of both types of mergers is similar in the range 
10"^ < Msub/Mhost < 10"^. Unfortunately, at this point 
our results from low mass haloes become limited by resolu- 
tion (i.e. we cannot identify smaller substructures) and the 
results from high mass haloes become dominated by Poisson 
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Figure 6. The number density of subhalo-subhalo mergers rel- 
ative to the mean density of subhaloes within r200 as a function 
of the distance to the centre of the host halo. The results from 
the MS are shown by solid lines while the results from the HS are 
shown by dashed lines. In each subpanel the dotted lines show the 
radial distribution of all subhaloes (regardless of whether they are 
merging or not) in the MS. Mergers involving subhaloes resolved 
with at least 50 particles are included in the plot. 



noise (i.e. less than one merger event in the whole simula- 
tion). Over the range that is reliably covered, we can see no 
strong systematic differences in Fig. \S\ between the resuhs 
derived from host haloes of different masses. This agree- 
ment is quite remarkable given the relatively large dynami- 
cal range resolved in the simulations. 

A merger between two objects is not always a straight- 
forward quantity to define in numerical simulations. The 
problem originates from the fact that any definition is in- 
trinsically linked to the mass and time resolution of the 
simulation. For instance, if in a higher resolution simula- 
tion we identify the remnant of a subhalo down to a smaller 
mass threshold, then the mass ratio of the merger, as well 
as the time at which it happens, could, in principle, disagree 
with the values measured in a lower resolution simulation. 
Similarly, with better time resolution, one could follow the 
mass loss of a subhalo more accurately which, in principle, 
could also change the measured mass ratio of the merger. 
To avoid these problems, we have chosen to use in Fig.[5]the 
mass of the satellite before accretion, rather than the mass 
at the moment of the merger. 

For all these reasons it is very important to note the 
agreement in Fig. [5] between the results from the MS (solid 
lines) and those from the HS (dashed lines). This agree- 
ment gives us confidence that our results are not sensitive 
to mass resolution. (Note that this is not true for subhaloes 
resolved with fewer particles as shown in the appendix.) Fur- 
thermore, the weak dependence of the quantities plotted in 
Fig. [5] on host halo mass confirms this conclusion. In prac- 
tice, a subhalo of Afsub/Mhost = 0.1 in a host of 10^^ h-^Me 



exhibits the same behaviour as a subhalo of the same frac- 
tional mass but in a halo of lO'^* h~^MQ even though the 
latter is resolved with 100 times more particles. This is quite 
remarkable. 

One of the reasons for the insensitivity to mass resolu- 
tion comes from our definition of a merger (see ij2.2p . We do 
not tag an event as a merger when we cannot identify the 
subhalo anymore, but rather when it has lost a fixed frac- 
tion of its most bound mass. This definition responds more 
to dynamical processes than to numerical ones. 

The implications of discrete time measurements are less 
clear for our definition of a merger. As an example, consider 
the case of very poor time resolution, and a halo that is 
just about to fall into a larger one. If tidal forces stripped 
off more than 95% of its mass before the next snapshot, 
then we would have identified this event as a merger. On 
the other hand, if the time resolution were good enough, 
we could have identified the subhalo at intermediate stages, 
updating its mass and the corresponding most bound 10 
percent. As long as stripping does not occur on a timescale 
much shorter than the time resolution, it is even possible to 
imagine that the line of descendants continues indefinitely. 
However, since a merger is not a discrete event, better time 
resolution does not necessarily imply a more accurate de- 
termination of a merger. With infinite time resolution, we 
would follow most of the merging process down to the point 
when mass resolution becomes important, i.e. every subhalo 
disruption would be caused by lack of mass resolution. 

However, the typical timescale for dynamical friction 
and tidal disrupti on is Tfric ~ tn for Maub/Mhost ~ 0.1 — 0.2 
I Jiang et al.ll200^ ). i.e. much longer than the time spacing of 
our simulation outputs (~ SOOMyr). Furthermore, subhalo 
mergers seem to take place very fast. Both these factors 
suggest that time resolution is not an issue for this study. In 
fact, we have checked that our results do not change if we 
choose snapshots that are twice as widely spaced as those 
used to build the merger trees. Nevertheless, we advise the 
reader to keep these limitations in mind. 



4.2 Characterization of subhalo-subhalo mergers 

In most cases the subhalo-central merger occurs very close 
to the potential minimum of the host halo. The spatial lo- 
cation of satellite-satellite mergers, on the other hand, has a 
very distinctive distribution. In the following subsection we 
investigate this further. 



4-2.1 Radial distribution of satellite-satellite mergers 

First, in Fig. |S] we look at the spherically averaged radial 
distribution of satellite-satellite mergers. The figure shows 
the number density of mergers, relative to the mean density 
of subhaloes within r2oo, as a function of the distance to 
the centre of the halo. We also display, as dotted lines, the 
distribution of all the substructures from the M^|. 



^ At first sight, the distribution of all substructures seems to dis- 
agree with the results of Section l3.3l Since the subhalo population 
is dominated by small mass objects, one would naively expect the 
distribution of all substructures to follow that of the smallest sub- 
haloes; as seen in Fig. |3] this has a slope which is always negative. 
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Figure 7. Probability distribution of the cosine of the separa- 
tion angle 6 between the progenitors of two substructures that 
are going to merge. The separation angle is measured at the last 
snapshot in which the subhaloes were identified outside the halo 
that hosts the satellite-satellite merger. Lines of different colours 
indicate mergers happening in haloes of different masses as indi- 
cated in the legend. 
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At every redshift plotted, the radial distribution of 
mergers is proportional to the radial distribution of sub- 
haloes. This implies that most of the mergers between sub- 
haloes do indeed occur in the outer regions of the host halo. 
Note that in these regions the background density is lower 
than in the inner regions, making it easier to identify sub- 
haloes. For this reason, we can follow satellite-satellite merg- 
ers down to structures resolved with 50 particles, as op- 
posed to the minimum of 200 particles we require for central- 
satellite mergers. 

Our results do not appear consistent with the naive ex- 
pectation from a gas kinetic theory approach that the num- 
ber density of mergers should be proportional to the number 
density of subhalo pairs, i.e Rmerge oc nlub - This discrepancy 
indicates that most of the satellite-satellite mergers do not 
occur because of random encounters between two unrelated 
substructures. We investigate this idea further in the fol- 
lowing subsection, where we look back at the orbits of the 
subhaloes that merge. 

4-2.2 Orbits of merging satellites 

Fig.[7]shows the distribution function of the separation angle 
9 between the progenitors of subhaloes involved in a merger. 
The angle is measured from the centre of the host halo in 
which the merger is going to take place, at the last snapshot 

However, in practice, the dominant effect is the high abundance 
of low mass host haloes in which only massive substructures can 
be resolved. As a result, the distribution of all subhaloes in the 
MS resembles the distribution of the most massive substructures. 



Figure 8. Three representative examples extracted from the 
MS for each of the two most common configurations between 
two satellite subhaloes that merge. The plots show the relative 
positions of the host and satellite halos in a time sequence, with 
time increasing from left to right. The black circles correspond to 
the halo that hosts the merger, while red and green circles show 
the positions of the satellites involved in the merger. The circles' 
radii are proportional to the half mass radius of each substructure. 
Class 1 : in this case the satellites were part of two separate haloes 
(red and green circles) before they were accreted into a larger halo 
(black circles). Class 2: both substructures belonged to the same 
halo before it was accreted into the larger structure which hosts 
the merger. 



in which both subhaloes were identified outside the halo that 
later hosts the satellite-satellite merger. It thus represents 
the angle between the subhaloes at the time they fall into 
the host halo. The first point to note is that the distribution 
seems to be universal in the sense that it is roughly indepen- 
dent of the mass of the host halo. (We have also checked that 
it is roughly independent of redshift.) However, the most im- 
portant feature is that the distribution is clearly dominated 
by small separation angles. About 65% of the mergers occur 
between subhaloes that were separated by less than 30 deg at 
the moment of accretion. (This percentage increases to 73% 
for an angle of 43 deg.) This demonstrates that the mergers 
are mostly between two or more systems that were already 
dynamically associated before they fell into the larger sys- 
tem. I f the gas kinetic theory approach of iMakino fc HutI 
(| 19971 ) applied to this case, then the mergers would be be- 
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tween subhaloes on random orbits, and we would expect 
a more uniform distribution in cos 6. (It would not be com- 
pletely uniform since the subhalo population is not isotropic, 
as shown in Fig. |4l) 

More information about the orbits of merging sub- 
haloes is given in Fig. [H] where we display three representa- 
tive examples of the two most common configurations of a 
satellite-satellite merger. These examples correspond to real 
sequences found in the MS. The plot tracks the position 
of substructures up to the snapshot of the merger (which 
happens at the rightmost position), starting on the left, 9 
snapshots earlier. We show as a black circle the halo that 
hosts the satellite-satellite merger and, as green and red cir- 
cles, the progenitors of the subhaloes involved in the merger. 
The red circle at the end of the sequence indicates the sub- 
halo resulting from the merger. The radii of the circles are 
proportional to the half-mass radius of the subhalo. 

The two most common configurations are as follows: 
Class 1: the progenitors of the subhaloes correspond to two 
separate haloes which were accreted at approximately the 
same time. Note that, as shown by Fig.[71 these haloes were 
spatially close at the time of accretion. Class 2: the merger 
occurs between two substructures that were part of the same 
halo before it fell into the host halo. In other words, there 
is a halo that contains two substructures which survived 
the accretion and subsequently merged. The merger event 
which started outside the main halo is completed inside it, 
as a subhalo-subhalo merger. 

Most subhalo mergers occur between substructures that 
are accreted close together both in time and location. Gen- 
erally, they were already part of the same system before it 
was accreted into a larger one, or were part of two separate 
haloes that were about to merge. This is probably a requisite 
for a subhalo merger to occur. The potential generated by 
the other satellite has to be at least comparable to that of 
the main halo. Hence, satellites accreted at difi'erent angles 
will follow relatively independent dynamical histories and 
are much less likely to merge. 
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Figure 9. Contour plot showing the logarithm of the number 
of satellite-satellite mergers as a function of the masses of the 
merging subhaloes at z = 0.5. The x-axis indicates the mass of the 
smaller subhalo while the y-axis indicates the mass of the larger 
subhalo. The different panels show the results for host haloes of 
different masses as indicated on each panel. The numbers in the 
bottom right show the number of mergers displayed in each panel. 
The vertical dashed lines indicate the 200 particle limit and the 
diagonal lines correspond to a 1:1 ratios between the masses of 
the two subhaloes. The horizontal lines show the mass limit on 
the more massive participant imposed by the choice of mass bin. 



4-2.3 The mass ratio of subhalo mergers 

In Fig. [51 we inspect the relative masses of the satellite sub- 
haloes which merge. The a;-axis indicates the mass of the 
smaller subhalo and the j/-axis shows the mass of the larger 
one. Interestingly, we find that, for the range of host halo 
masses plotted, the most common merger is that between 
two substructures of dissimilar masses, Msub,i ~ 10xMsub,2. 
Note that this trend is contrary to the naive expectation 
whereby the mergers are simply proportional to the abun- 
dance of substructures, in which case the maxima would 
be located around the line Afsub,i = Msuh,2- However, it 
is roughly consistent with the idea that substructure merg- 
ers happen between two structures that were part of the 
same halo before accretion. For instance, if the most com- 
mon merger happens between the main subhalo and its most 
massive substructure, then, as we have seen, we would ex- 
pect to find a mass ratio of 1:25 (see Fig. [2]) and the maxima 
of Fig. [9] along Msub.i ~ 10 x Msub,2 



4.3 Merger probability since accretion 

Finally, in Fig. [10] we plot the fraction of subhaloes at a 
given redshift that have had a merger with another satellite 
subhalo since accretion into the current host halo. The top 
panels show mergers between satellites with a mass ratio 
greater than 0.03, i.e. in which the less massive subhalo has, 
at least, 3% of the mass of the larger one. In the bottom 
panels we consider mergers between subhaloes with more 
similar masses: the minimum mass ratio is 0.3. 

The fraction of current subhaloes which have experi- 
enced a merger in the past is a quantity strongly affected by 
resolution. For instance, in the history of a subhalo resolved 
with 1000 particles, because of our 200 particle mass cut on 
subhaloes, we can only record mergers with other subhaloes 
which account for at least one fifth of the final subhalo mass. 
On the other hand, if our current subhalo is resolved with 
10000 particles, then a much wider range of merger mass 
ratios can be tracked. These considerations are further com- 
plicated by the fact that we expect the measured mass of 
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Figure 10. The fraction of substructures that have experienced a merger with another substructure since the time of accretion into the 
current host halo. The x-axis gives the subhalo mass at the redshift shown, while the ratio M2/M\ on the y-axis is for the two progenitors 
of the subhalo at the time they merged. The results from the MS arc shown by solid lines while the results from the HS are shown by 
dashed lines. The coloured lines represent the data from haloes of different masses, as indicated by the key. The two rows correspond to 
different mass ratios between the subhalo progenitors involved in the merger: A/gub 2 > 0.03 Mgut 1 (top row) and Msub,2 > O.SAfgui, 1 
(bottom row) where Msub.2 refers to the larger satellite involved in the merger. The three panels display the results for substructures 
identified at redshifts z = 1, 0.5 and respectively. 



a subhalo to be less than the mass of its progenitors at in- 
fall, due to tidal disruption and stripping; hence an object 
that is below our 200 particle limit at a particular redshift 
could have been above this mass cut when it experienced 
the subhalo-subhalo merger. 



To improve statistics, whilst at the same time attempt- 
ing to avoid building a resolution dependence into our re- 
sults, we relax the particle number constraint on subhaloes 
for this exercise. At the redshift a subhalo is identified (i.e. 
the redshift plotted in Fig. [10}, we consider subhaloes of 
30 particles or more. At the redshift of the subhalo merger, 
the progenitors must both have 50 particles or more to be 
counted. 



Fig. [TD]shows that the probability of a subhalo merger is 
constant for subhaloes of different mass. About 1% percent 
of subhaloes have had a merger with another subhalo with 
a mass ratio > 0.3. For a mass ratio > 0.03, this fraction 
increases to ^ 10%. We also note that these fractions show 
a weak decrease with redshift. 



5 SUMMARY AND CONCLUSIONS 

We have used the Millennium simulation, together with a 
simulation which has 10 times better resolution but about 
100 times smaller volume, to investigate the general proper- 
ties of the substructures within dark matter haloes, includ- 
ing their merger rates. Our main findings can be summarized 
as follows: 

In agreement with previous studies, we find that the 
mass function of low and intermediate mass subhaloes fol- 
lows roughly a power-law. However, we also find an expo- 
nential cut-off in the mass function at high subhalo masses. 
We have provided an expression, Eq. [2j that describes this 
behaviour accurately. We also detect a small but system- 
atic dependence of the number of subhaloes on the mass 
of the host halo. On average, at a given fractional mass, 
Afsub/Mhost, high mass haloes contain more low and inter- 
mediate mass substructures than their less massive counter- 
parts. In contrast, we find evidence that high mass haloes 
contain fewer high mass subhaloes than do low mass haloes. 
In spite of this, the fractional mass of the first, second and 
third most massive substructures is insensitive to the mass 
of the host halo and of the redshift. 

We confirm that the radial and angular distributions 
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of subhaloes are roughly independent of the host halo mass 
and redshift. However, we find that the radial distribution 
does depend on the subhalo mass relative to that of the host 
halo. The subhalo distribution is less concentrated than the 
dark matter, but the radial distribution of low mass sub- 
haloes tends to be more concentrated than that of high 
mass subhaloes. This difference can be understood as re- 
sulting from the different efficiency of dynamical friction in 
subhaloes of different mass. On the other hand, these dis- 
crepancies between the radial distributions of low and high 
mass subhaloes disappear in the outer parts of the halo, as 
seen i n recent ultra-high r esolution simulations of galactic 
halos (jSpringel et al.ll2008l ). 

The angular distribution of subhaloes tends to be 
aligned perpendicular to the spin axis of the host halo. This 
is probably due to an anisotropic mass accretion - merg- 
ers happen preferentially along filaments. The alignment is 
strong for the most massive subhaloes, but is much weaker 
for low mass substructures since, on average, they have spent 
a few orbital times inside the halo which would randomize 
their orientation. 

We have found that satellite-satellite mergers do occur. 
Over most of the mass range resolved in our simulations, 
they are subdominant when compared with mergers between 
satellites and the central subhalo. However, we see some 
indication that satellite-satellite mergers are equally likely 
to satellite-central mergers for the lowest mass subhaloes 
(Msub/Mhost < 10~^). As for many other subhalo properties, 
the merger rates appear to be a function of the fractional 
subhalo mass only, and are independent of the particular 
host or subhalo mass. 

The radial distribution of satellite-satellite subhalo 
mergers closely follows the radial distribution of subhaloes. 
This implies that most of the subhalo mergers happen in the 
outer layers of the halo. For the most part, these mergers 
involve subhaloes that are already dynamically associated 
before accretion into the main halo, i.e. they were either 
part of the same halo, or of two separate haloes that were 
accreted at similar times and locations. At every redshift, 
most of these subhaloes which subsequently merged were 
closer together than 30 deg as seen from the centre of the 
halo that hosts the merger, at the time they fell in. 

Finally, we find that a small fraction of the high-mass 
subhaloes has experienced a merger with another subhalo 
since accretion into the current host halo. The values depend 
on the mass ratio of the merger, but vary from a few percent 
for mass ratios greater than 0.3 to ~ 10% for mass ratios 
greater than 0.03. 

In spite of using some of the largest simulations to date, 
our results could still be affected to some extent by numeri- 
cal resolution. Due to the rarity of the events we are trying to 
study, it is difficult to find a range of substructure and host 
halo masses where we have, at the same time, (i) enough par- 
ticles to resolve substructures well, (ii) enough haloes to dis- 
tinguish real trends from cosmic variance, and (iii) enough 
subhaloes to establish their properties and dynamics. Fortu- 
nately, as we have shown, many properties can be described 
as a function of only the fractional subhalo mass. In these 
cases we are observing the same system resolved with many 
different numbers of particles, so it is reassuring that we find 
the same trends for different host halo masses. This gives us 
confidence that these results are robust. On the other hand. 



quantities which scale with halo mass are much less reliable 
and could still be affected by resolution effects. Much larger 
simulations, currently beyond reach, will be needed to check 
them. 
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APPENDIX: NUMERICAL EFFECTS 

Numerical artifacts can pose serious problems in obtaining 
a robust estimate of various properties of the population of 
subhaloes. For instance, two-body encounters, particle heat- 
ing, or force softening could easily dilut e substructures tha t 
are not resolved with enough particles l|Moore et al.|[l99^ . 
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Figure 11. The mean number of satellite-central subhalo merg- 
ers per subhalo and per unit of time as a function of the subhalo 
mass. The solid lines show the results from the MS while the 
dashed lines show the result from the HS. The coloured lines rep- 
resent the results from haloes of different mass, as indicated by 
the legend. Note that we display results from subhaloes with 20 
particles or more. The upturn in Nm for low mass subhaloes is 
due to the inclusion of subhaloes resolved with fewer than 200 
particles. Once the > 200 criterion is applied, the upturn dis- 
appears as shown in Fig [5] 



These problems translate into an overestimation of the num- 
ber of satellite-central subhalo mergers in each timestep. 

Such a feature is clear in Fig. 1111 which is similar to 
Fig. O but for satellite-central mergers only and including 
subhaloes with less than 200 particles. For these objects, 
we can see a strong disagreement between the merger rate 
of substructures in the simulations with different resolution 
which is manifest as an upturn in the curves. However, the 
upturn disappears for subhaloes with A'' > 200 which is the 
limit set in this paper. 

An overestimation of the destruction rate also has impli- 
cations for other quantities such as the abundance and radial 
distribution of subhaloes. For instance, the subhalo mass 
function shows a cut-off at low masses compared with the 
expected power-law behaviour when we include subhaloes 
resolved with fewer than ~ 50 particles. (This quantity is 
less affected since most of the subhaloes are in the outer 
layers of the halo.) On the other hand, the inner part of the 
radial distribution is more sensitive to these effects. Once 
subhaloes with fewer than 200 particles are included in Fig. [3] 
the distribution becomes less centrally concentrated. 

Our convergence study indicates that 200 particles is 
the limit below which results are unduly affected by reso- 
lution. This is why we have adopted this minimum particle 
count throughout this chapter, except when otherwise stated 
explicitly. This choice should minimize finite-resolution ef- 
fects. 



